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The existence of a spin disordered ground state for the frustrated Ji — J2 — J3 Quantum Heisenberg 
Antiferromagnet on a square lattice is reconsidered. It is argued that there is a unique action which 
is continuous through the whole phase diagram, except at the Lifshitz point, so that the Neel and 
helicoidal states can not coexist and that there has to be an intermediate spin liquid state. To 
show it, a detailed study combining Spin- Wave Theory, Schwinger Bosons Mean Field Theory and 
Oh. a scaling analysis of the appropriate hydrodynamic action is performed. When done carefully, all 

these theories agree and strongly support the existence of the spin liquid state. 
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I. INTRODUCTION 
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The past few years have seen a considerable effort in 
understanding the nature of the ground state and exci- 
tations of low-dimensional Quantum Antiferromagnets, 
due to its close link to high-T c superconductors jl]|],|3],|| ■ 
One of the most fundamental discussed issues is the ex- 
istence and properties of spin liquid states. 

The way in which one can induce a spin disordered 
state in a magnet is by increasing fluctuations about the 
O ■ ordered state, for instance, by adding frustration to the 
bare Heisenberg model -antiferromagnetic coupling con- 
stants to next to nearest neighbors- [|| . We study in this 
paper the phase diagram of the quantum J\ — J2 — J3 
model at zero temperature (all Ji > 0). We show in fig- 
ure 1 the phase diagram of the corresponding classical 
model. §. We see that there are four possible states 
-Neel, collinear and two helicoidal- separated by bound- 
aries. The transition between Neel and collinear states is 
discontinuous because they have different symmetry. The 
same happens in the boundary between the (<?, q) and 
(q, 7r) helicoidal states. The other transitions (marked 
with dashed lines) are continuous and separate states of 
the same symmetry. The Classical Critical Line (CCL) 
J\ — 2J2 — 4J3 = is the separatrix between Neel and 
helicoidal states, and its endpoint, J3 — , is a Lifshitz 
point. 

The first studies on the phase diagram of the quantum 
model were performed by Chandra and Doucot lq] (for 
J3 = ) using Linear Spin Wave Theory (LSWT) and 
by Ioffe and Larkin |lC|], who made a Renormalization 
Group (RG) analysis of a generalized nonlinear a model. 
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Both of them concluded that the combined effect of quan- 
tum fluctuations and frustration were strong enough to 
destroy long-range-order (LRO) along some critical line 
yielding a spin liquid state. Later on, Moreo et al || 
and Chubukov generalized the LSWT analysis to 
J3 7^ (see figure 1). This scenario was supported by 



(for J3 = ) and a momentum-shell RG study |15||. On 
the other hand, several mean field-BCS-type theories of 
interacting bosons have challenged this result, either for 
J 3 = 0,|l^,|l|,|l9| or for the general case @. They 
predict that the stability of the Neel state is enhanced 
over the classical case due to quantum fluctuations (Vil- 
lain's order from disorder pl[) along the whole critical 
line, while that of the spiral state is not modified in an 
essential way. 

The purpose of this paper is to clarify this controversy. 
Our main points will be: 

1. The CCL separates states with the same symme- 
try. The transition is continuous classically, in the 
sense that the parameters of the long wavelength 
action S[Ji] vary smoothly along the transition. 
We argue that the same happens for the quantum 
case, S[Ji, S] . In particular, we show that, although 
there is an enhancement of the stability of the Neel 
state close to the CCL, there is a bigger decrease 
of the stability of the spiral state. This means that 
the CCL is continued in a Quantum Critical Plane, 
QCP, which depends on S and is tilted upwards 
with respect to the plane J\ — 2 J 2 — 4J3 (see figure 



2. Because the quantum action is continuous through 
the QCP, we can use a generalized Sigma model 
to describe the physics of the transition along the 
QCP and at its Neel side. We evaluate that ac- 
tion along the QCP using a Renormalization Group 
analysis and find that the relevant coupling con- 
stant always flows to strong coupling meaning that 
the spin state is disordered. We then proceed to 
show that the spin liquid state exists not only along 
that plane but also in a finite region around it. 
As a corollary, we link this scaling analysis for the 
strongly frustrated magnet with the one performed 
in reference 2, for the weakly frustrated case. 

3. The end point of the CCL is a Lifshitz point. In this 
case, the counterpart of the Neel state is not spi- 
ral, but collinear. Because both states have differ- 
ent symmetry, the quantum action for both states 
is different, and therefore we cannot establish any 
continuity principle. We discuss this case in the 
last section. 
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II. MICROSCOPIC STUDY OF THE Ji - J 2 - J 3 
MODEL 



The purpose of this section is to study the stability 
diagram of the quantum J\ — J2 — J3 model around the 
CCL. Because the whole line has the same physics, we can 
concentrate our analysis just to one, representative point 
of the phase diagram. We will choose the plane J2 = 
for simplicity, in which case the CCL is J3 = 0.25Ji. We 
will work in units in which J\ — 1. 

The energy of the classical J\ — J2 — J3 model is given 

by 

- - S 2 

E = X] Ji i Si ' S i = "jf Jq °' c! ^ 
<i,j> 

where Jq is the Fourier transform of «7y and Qo.ci is the 
pitch wave vector of the spiral state, which is obtained 
by minimizing the energy: 8qJq cl = 0. Qo.cZ is a con- 
tinuous function of J2/J1 and J3/J1, with a jump in its 
derivative at the critical value J\ — 2 J 2 — 4 J 3 = 0. This 
equation defines the Classical Critical Line. The classical 
spin stiffness is given by 

, _ d 2 E 

and it also varies continually with Qo,d- 

We can then write the most general action for an heli- 
cal magnet by noting that the order parameter space is 
0(3) x U{l)/U{\) llHEl 

S = ~\j d 2 xTr{A l P l A l } = ~\ J d 2 xp a (A«) 2 (3) 

where A, = A°; T a = g~ 1 dig is a pure gauge field, 
g(x)eSO(3) and T a e Lie[50(3)]. The gauge field is 
equivalent to a twist of the order parameter and serves to 
define the spin stiffnesses, p a = (p x , p x , p z ). The action 
S[p(Qo,rt)] is unique for both Neel and helical states and 
the spin stiffnesses are continuous throughout the whole 
phase diagram. 

We argue that, after introducing quantum fluctuations, 
there is still a unique action for the whole phase diagram, 
which varies continually when passing from the Neel to 
the spiral state through the Quantum Critical Plane. The 
action picks up an extra piece due to the fact that the 
fields become dynamical and an implicit dependence on 
S through Qo{S): 

S = -\J dx d 2 x( Xa (A a ) 2 + Pa (A«) 2 ) (4) 

All this implies that if the stability of the Neel state 
gets enhanced due to quantum fluctuations, (1) the pitch 
wave vector of the helical state also renormalizes so as to 
adjust smoothly to the new boundary of the Neel state, 
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where its value is (tt, 7r), (2) the spin stiffness of the heli- 
cal phase is equally decreased, so that both go to zero at 
the same point. To substantiate this ideas, we will per- 
form in this section a Spin Wave Theory (SWT) study 
of the magnetization, the pitch wave vector and the spin 
stiffness and complement it, eventually with results from 
Schwinger Bosons Mean Field Theory (SBMFT). Details 
of the derivation and notation are relegated to appendices 
A and B. 

As a starting point, let's determine the phase diagram 
from LSWT and SBMFT. The magnetization for the he- 
licoidal Qq — ((fa, So) and Neel Qo = (tt, it) states from 
LSWT is 

<C*> Q I \ 1 y J fc + J ± ~ 2 J Qo f5l 

<b >— b+ 2 2N 2_, E{k) (5) 

Qo is obtained by minimizing the ground state en- 
ergy: cos(go) ~ cos(<7o, c/) = —477 (f° r the Neel state, 
cos(<5o) = —1). The denominator is the energy of the 
spin waves. The value of S for which the staggered mag- 
netization < S z > is zero determines the boundaries of 
the long range ordered states. 

We plot in figure 2 the stability diagram we have ob- 
tained using LSWT and SBMFT. The theories are in dis- 
agreement: (l)LSWT predicts a decrease of the stability 
of both ordered states, with a line of second order phase 
transitions ending at the critical point; (2) SBMFT pre- 
dicts an enhancement of the stability of the Neel state. 
This enhancement is also seen in numerical diagonaliza- 
tion studies |jj . On the other hand, the boundary of the 
spiral state doesn't move towards the right as we intro- 
duce quantum fluctuations. This is a spurious result, due 
to the fact that large-N theories don't give the correct 
Goldstone mode structure for the spin wave spectrum. 
As a result, we obtain a first order phase transition. 

Let's concentrate for a moment on the Neel state and 
compute the expression for the staggered magnetization 
close to the classical frustrated point. Evaluating the 
leading divergences we obtain 

<S»>~5+i-oln(pcj) + |— + 0(4) ( 6 ) 
2 b pd b z 

where p c i = J\— 4 J3 is the classical spin stiffness which we 
use here as a cutoff for the infrared divergent integrals. 
This series can be resumed: computing £ (k) -the spin 
wave energy to next to LSWT order in 1/S-, we find: 

<5»> ~S+i- 7 ]n(p(S)) (7) 

p(S) ~ pd - -^(0:3 -aa) 

which is the LSWT result with a renormalized spin stiff- 
ness. It serves us to define the QCP, i.e.: the line where 
the spin stiffness, computed to next to LSWT order is 
zero. As shown in figure 2, p(S) is tilted to the right of 
the classical value, with a slope ~ 5.2. This fact supports 
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the order from disorder conjecture from the point of view 
of SWT. 

Based in our symmetry argument, we argue that this 
line has to be the same for both the helicoidal and the 
Neel states (at least for large S). In other words, the 
enhancement of the stability of the Neel state is accom- 
panied by a similar reduction of the stability of the heli- 
coidal one so that both states match each other continu- 
ously. To prove it, we compute the LSWT correction to 
the pitch wave vector and to the stiffness: 



% %,cl 



2^ (2) j Qo , c! V \J±-Jqo,J 



s2 dj J Qo,.i 
2 

S 



+ 2 { d< ? )j Q a ,ci ~ (sin(©,ci) + 8J3sin(2gb,ci))5A?5 

1 V (2Jq , e , -Jk- J±)dl 2) J Qo , cl - (J Qo . cl - W (2) J± ' 
2^ E(k) 



•0(1) 



The spin stiffness of the unfrustrated magnet is p = 
Pci(l-^r) (= 0.764p c! for spin 1/2). This estimate is in 
accordance with that obtained using the renormalization 
of the spin wave velocity and the susceptibility H, p = 
^,„(l-f)(l + 2xf). 

We draw in figure 3(a) the LSWT correction to q c i for 
the spiral state to show that it matches with (n, 7r) con- 
tinually, although the boundary of stability of the Neel 
state be shifted to the right. Figure 3(b) is a plot of the 
LSWT correction to the spin stiffness, and shows that it 
is enhanced for the Neel state close to the critical point, 
signaling the order from disorder phenomenon. The im- 
portant point is that it is decreased by a bigger amount 
on the other side. We stress that this beautiful results 
are already obtained in LSWT. 

Before entering in the next section, we would like to 
summarize the most important conclusions we have ob- 
tained so far. We can very safely state that the stability 
of the Neel state is enhanced over the classical values over 
the whole critical line, and that this is due to the fact that 
the spin stiffness p(S) gets renormalized upwards due to 
quantum fluctuations, so that it defines a QCP tilted to 
the right of the CCL. The stability of the helicoidal state 
is reduced in the same way, so that both states match 
smoothly and the long wavelength physics is continuous 
through the transition point. Because of that, a modi- 
fied a model which describes the spin wave excitations of 
the strongly frustrated Neel state is appropriate to study 
the QCP and its neighborhood. These conclusions don't 
apply to the J\ — J2 model due to the fact that in that 
case the two states have different symmetry and the spin 
stiffnesses close to the tricritical point are different. 
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III. SCALING ANALYSIS 



Wc devote this last section to perform a Renormal- 
ization Group study of the hydrodynamic action appro- 
priate for the strongly frustrated magnet. It is a gener- 
alization of the 0(3)/C/(l) = S 2 quantum non linear a 
model. We will review in section 4.1 well-known results 
for the conventional 0(3) case, published by other au- 
thors, which are pertinent to our discussion. In section 
4.2, we pass to the strongly frustrated case, and show that 
along the QCP the system always flows towards strong 
coupling. The results of this section complement our pre- 
vious findings from microscopic models and give us a very 
consistent picture of the phase diagram of the J1 — J2 — J3 
model. 



A. Unfrustrated case 



It is generally believed that the long-wavelength ac- 
tion of a non-frustrated Quantum Antiferromagnet in d 
dimensions is given by a (d + l)-dimensional non- linear 
a model [f25| 



S=£ [°° dt [ d d r {{Vnf - \(d t n) 2 

^ J-00 Ja I c 

dt [ d d f{{Vh) 2 - (din) 2 } 



(9) 



1 
2^o" 



where f = r/a, t — tc/a are dimensionless variables, 
Po = JiS 2 a 2 ~ d is the bare spin stiffness, Co = 2^/dJiSa 
is the spin wave velocity and 
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1 _ 2\fd 
Poa^ 1 V/W^ 1 s 



is the bare coupling constant. All these magnitudes are 
defined at the scale of the lattice cutoff, a. 

In 1+1 dimensions, one has to add to this action a 
topological 6 term, which distinguishes between integer 
and half-integer spins p5|] . It has been proven that such a 
term gives a zero contribution in 2+1 dimensions for any 
topological configuration of the n-fields |2(||27],^8| . More- 
over, as long as the n-field is continuous at all points 
no Hopf term appears in the action. Once this condi- 
tion is relaxed, though, that term appears and affects 



the physics of the system 26 29] . Such a situation could 



arise if the a model were in its disordered phase. 

One of the main issues concerning the nonlinear quan- 
tum a model is the connection with the Classical one, 
which is much better known. This connection can be 
done in principle when the contribution of the Berry 
phase is zero. Strictly speaking, the mapping is rigorous 
only for large S and long wavelengths. We therefore face 
two problems when we try to link it with the microscopic 
quantum models of interest: 
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• We don't know by certain if the action is still valid 
in the extreme quantum limit. In this respect, nev- 
ertheless, we can argue that it gives the correct spin 
wave spectrum and interaction between spin waves 
at long wavelengths, as one can check by calculating 
the staggered magnetization, spin wave energy and 
damping,. . . in the first few orders of perturbation 
theory 

• A more difficult point is the identification of the pa- 
rameters entering in the action in terms of the ones 
appearing in the microscopic model. The prob- 
lem is that the values obtained in the derivation 
are well suited only for large spin and long wave- 
length scales, and they don't have to hold necessar- 
ily for small S. In fact, they don't. One possibility 
is to connect those values with the ones obtained 
by SWT. But it is not clear which order in 1/5 
should one use as input of the continuum model. 
That should depend on how strongly the Goldstone 
modes interact and, therefore, on the dimension. 

Polyakov pioneered the work on the Classical NLcrM 
by recognizing that close to its lower critical dimension 
Goldstone modes interact strongly and become critical. 
He proceed then to solve the problem using a Renormal- 
ization Group approach. Later on, the same problem was 
addressed by Brezin and Zinn-Justin |3l}| in the context 
of the e-expansion (e = d — 2) to order e 2 and by Nel- 
son and Pelcovits |32(| using momentum-shell recursion 
relations Q to one loop order, i.e.: to first order in e. 
The approximation is expected to be very good in two 
dimensions, but it is not so clear that this is the case in 
three (where e = 1). 

Once we have an action to describe the system in 
2+1, we have to distinguish between short wavelength 
quantum critical fluctuations and the long wavelength 
antiferromagnetic magnons. The crossover length be- 
tween both is the Josephson correlation length £j 
£j tells us the extent to which the quantum model can be 
mapped to the classical one. This problem was addressed 
by Chakravarty, Halperin and Nelson in an important 
publication . CHN combined the techniques of Nelson 
and Pelcovits with the method developed by Hertz pq ] 
to deal with quantum critical phenomena at finite tem- 
peratures in two spatial dimensions. They showed that: 

• At T — the model is equivalent to a 3-dimensional 
classical one where g ~ 1/5 plays the role of the 
temperature. In that case, the QNLerM has a phase 
transition from a Neel ordered to a (quantum) dis- 
ordered state, at a value g c of the coupling constant. 
The Josephson correlation length turns out to be of 
the order of a few lattice spacings for g < g c , but 
grows with g until it diverges at g c . 

• For g > g c the system is in a quantum disordered 
phase. In that phase, CHN obtain £ ~ s^-s ^ or ^ ne 
correlation length. It is very likely that instanton 
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contributions to the Hopf term will stabilize spin- 
Peierls order here |2{J . 

• At finite, but not too large, temperatures in the 
region where the ground state is Neel ordered {g < 
g c ) we can integrate all quantum fluctuations to 
map the QNLcrM into a classical problem with a 
renormalized spin stiffness and a short wavelength 
cutoff of order Tic/kBT. 

Our goal is to analyze the stability of the Neel state 
when introducing frustration. We will therefore use CHN 
approach in the zero temperature formalism. To begin 
with, we recalculate some of the magnitudes already dis- 
cussed by the former authors. We defer to appendix C 
the derivation of the scaling equations, and quote here 
only the results. We obtain: 



(10) 



dg 






dl 


= -§ + 


f 


dp 






~dl 


= -pa 




da 






~dl 


= <?§ 






S c 


0.225 


9o 


s 


S 



where a is the component of n along the quantization 
axis and g = Integrating 

floe"' 

gi = z — r (11) 

1 -50(1 -e l ) 

Pi = po(l-<?o(l - e - *)) 

co = <7oo(l — go) 

We plot the flow diagram in figure 4(a). There are 
two different phases, separated by the straight line g = 
1. To the left of this separatrix (g < 1), we have the 
renormalized classical phase, where there is Neel long 
range order renormalized by quantum fluctuations. All 
the flow lines end at the p axis, which is a line of fixed 
points. The solution of the scaling equations is: 

9oo - (12) 

S 

Poo = po(i - go) = po(i - -j) 

§ 

< a > = l- 0o = l- — S <cr > = S - S c 

At the right of the separatrix, we have the quantum 
disordered phase, where interaction between magnons is 
so strong that confine them, generating an internal length 
scale I* — ln( 90 . ) in terms of which we define the cor- 
relation length 
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The coupling constant g diverges and the spin stiffness 
goes to at that length scale. 

Finally, we think it is worth remarking that although 
figure 4(a) shows flow lines for p ~ 0, that result is spu- 
rious, because in that case the action for the QNLcrM is 
empty and one has to take into account quartic terms in 
the gradient expansion of the Heisenberg Hamiltonian. 



B. J1-J2-J3 Model 



We are now in a position to understand better the 
strongly frustrated case. We will use a modified a model 
appropriate for both the weakly and strongly frustrated 
magnets in the Neel state. We will try to discern if the 
coupling constant flows towards strong coupling along the 
QCP, signaling that the magnet is disordered for what- 
ever the value of the spin. 

Close to the QCP the bare spin stiffness, p(S), is very 
small and we are forced to take into account effects due 
to quartic terms, just as we did in our calculation of the 
staggered magnetization in SWT. The action is 



S = - / dtd 2 r { X °(d t n) 2 - p°(S)((d x h) 2 ) + {dyhf)) (14) 

z J — oo J a 



a°a 2 

+ -^-{{d xx n) 2 + {d yy h) 2 ) - a%a 2 (d xx n)(d yy h) 



We take p°(S) as given by SWT (i.e.: the stiffness 
renormalized by "short" wavelength quantum fluctua- 
tions), instead of p — {J\ — 2J 2 — 4J 3 ) 5 2 , and use the 
values obtained in the gradient expansion for the rest of 
the parameters entering in the action: 

a? ~ ( ,h - 2J 2 - 16J 3 )5 2 



9° = 



8J1 



^pJS)^a V Ji ~ 2 J 2 - 4 J 3 S 



We can already see one of the problems which arise in 
the frustrated case: the coupling constant depends not 
only on S (as in the unfrustrated case), but also on the 
exchange interactions Ji, J2 and J3. Therefore, we don't 
have such a clear indicator about the magnetic order of 
the microscopic model as before. Moreover, go diverges at 
the QCP, and therefore is not well-suited for our analysis. 
This means that we will have to cook a charge which 
doesn't contain p(S) in the denominator. 

After performing the perturbative RG expansion we 
arrive at the following scaling equations: 

IT —»<» + <") 



9 



dc _ _ {<J\ - 2cr 2 ) 
dl ~ ° 9 8p 

where I is the loop integral 

24p - 2(7i + (ci + 6o- 2 ) sin 2 

This integral diverges at the Lifshitz point, signaling 
that the scaling analysis breaks down there. We see that 
this point is quite pathological, and has to be treated in a 
different way from the rest of the phase diagram. We have 
written down the scaling equation for c to show that in 
the frustrated case, it does get renormalized downwards. 
These scaling equations reduce to those obtained by CHN 
for the unfrustrated case and to those of Ioffe and Larkin 
when p(S) = 0. 

We will pass on now to the task of finding an unam- 
biguous indicator of the ordering of the system. We have 
two possibilities to obtain a charge which roughly go like 
1/5 and which give us the separatrix between the ordered 
and disordered states: 

1. We can choose Go ~ This kind of charges di- 
verge if po(S) is close to the QCP. They have a typ- 
ical behavior Q ~ — a(p) Q + j3{p) G 2 , which means 
that there is a separatrix but, because Go >> 1, 
one is always in the disordered phase. 

2. Go ~ y5=5j- I n this case, the system always flows 
towards strong coupling and there is no separatrix. 

There is a third possibility: we can evaluate the action 
strictly along the QCP and define, following Ioffe and 
Larkin |[(|, the following charge: 



We then assume that A = (<7ia 2 ) 1/4 I / p 1/2 (S) doesn't 
get renormalized. We have checked that the QCP is sta- 
ble under iteration of the scaling equations. We find that 
G satisfies: 

^ =C?2 :Gl = l-cTil-lo) (18) 

It flows towards strong coupling, which means that it 

stays effectively at its lower critical dimension, with the 

i „ 

generation of a correlation length £ ~ ae 7 ^ ~ ae . 

We should draw now the separatrix between the or- 
dered state and the spin liquid. Ioffe and Larkin argue 
that for finite p(S) the relevant action is not given by 
equation 14, but by 

i poo pL 

S = - / dtd\ { X \d t nf (19) 

a°a 2 

+ -L-((d xx h) 2 + {dy V fi) 2 ) - a°a 2 (d xx n)(dyyh) 
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where ln(^y) acts as a cutoff in the scaling equations. 
In that case, the separatrix is given by 

l + GoHp(S)) =0 (20) 

We stress that Ioffe and Larkin's action is different 
from the naive hydrodynamic one, equation 14. It is not 
clear to us which one of them really describes the physics 
of the Ji — J2 — J3 model close to the QCP. In any case, 
there is no doubt that the separatrix has to exist, be it 
an exponential or not. We show in figure 4(b) the flow 
diagram of the complete J\ — J2 — J3 model, where we 
draw a tentative separatrix for small p and join it with 
the one given for the weakly frustrated case. It is amusing 
to note that turning the figure by 90°, the diagram is the 
same as the one given by SWT and SBMFT, if we take 
into account that we are writing the dressed p(S) in the 
ordinate axis. 



C. Discussion 

We recapitulate here the main conclusions we have 
reached: 

• We have found that the system is at its lower crit- 
ical dimension along the QCP: the charge flows to 
strong coupling, signaling a spin liquid state. 

• There is a separatrix between the spin liquid phase 
and the Neel ordered state. That line begins at the 
point (p(S),g) = (0,0) and links the CHN separa- 
trix for sufficient large values of the spin stiffness. 

• The flow diagram looks very much like the stability 
diagram that we obtained from SWT and SBMFT, 
with an enhancement of the Neel phase due to order 
from disorder. 

• There is a dimensional crossover from the region 
p(S) << 1, where the system is 2+1 dimensional, 
and g' = —g + g 2 to the region where the spin stiff- 
ness is very small, and the system is at its critical 
dimension, d c , so that G = G 2 . 

IV. Ji - J 2 MODEL 

We devote this section to study the J\ — J2 model. We 
have already argued that it is the most complicated case 
of the whole phase diagram, and probably, the least rel- 
evant. We nevertheless study it because it has attracted 
a lot of attention. 

We begin, as in section II, by plotting the phase di- 
agram obtained from LSWT and SBMFT (figure 5(a)). 
The expressions for the magnetization of the Neel and 
collinear states from LSWT as a function of J2/J1 are 

§: 
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<s*> =s + l-ly - i-Mi-Tk) {Ned) (21) 

= S + l-iy - l + Acos^) 

2 JV ^ ^(1 + Acos(fc a )) 2 - cos^) 2 ^ + cos(fc a )) 2 



where A = 1/2 J 2 . The denominators are the velocities 
of the spin waves. The value of S for which the magneti- 
zation < S z > goes to zero defines the boundary we are 
looking for. 

As already stated by other authors , LSWT pre- 

dicts a finite region where the ground state is disordered 
around the classical tricritical point J<i = 0.5 for any S 
(and in particular, for 5 = On the contrary, mean 
field theories predict that the stability of the Neel state 
is enhanced by quantum fluctuations with a reentrant 
boundary joining the classical frustrated point with zero 
slope. This fact means that there is, ineludibly, a region 
of coexistence of Neel and collinear states and therefore, 
a finite hysteresis and a first order phase transition. This 
first order phase transition extends even for 5 < 1/2, so 
that there isn't a spin liquid state for spin one-half. 

We would like to devote some words to some interesting 
features of the SBMFT boundary of the collinear state. 
In fig. 5(b) we show the SBMFT results in a logarithmic 
scale. It shows a bump where the stability of the collinear 
state is even bigger than that of the unfrustrated case and 
then tends to the correct result for J2 — > 00 (1/5 ~ 5.1). 
In the same figure we show the results obtained in nu- 
merical diagonalizations ]lj| for spin one-half. We see 
that the boundaries of the ordered states from mean field 
theory are much higher than the results of numerical di- 
agonalizations. These two pieces of evidence support our 
belief that mean field theories overestimate the stability 
of the ordered states and that there is a spin liquid in a 
very small region at spin 1/2 for J2 ~ 0.6. 

We solve now the discrepancy between LSWT and 
SBMFT. We evaluate the staggered magnetization of the 
Neel state in the next to LSWT order in 1/5 close to the 
CCL: 

<^>^5 + i-aln 2 (p cJ ) + |^^ + 0(^) (22) 
2 5 pa S z 

where we use p c i = J± — 2 J2 as a cutoff for the infrared 
divergent integrals and a and (3 are constants. This series 
can be resumed by computing £ (k), the spin wave energy, 
to next to LSWT order in 1/5. Using the spin wave 
velocity 

P{S) ~ pel ~ ~cf( a 2 ~ "l) 

as a cutoff, instead of the spin stiffness, we find: 

<5*>~5+i- 7 m 2 (p(5)) (23) 

which means that 1/5 corrections to LSWT just renor- 
malize the CCL where the spin stiffness (and spin wave 
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velocity) is zero. The QCP for the Neel state along which 
the spin stiffness is zero gets tilted to the right as one goes 
further in perturbation theory in 1/5. Therefore, SWT 
also predicts an enhancement of the stability of the Neel 
state and therefore, order from disorder. In other words: 
it tends to coincide with SBMFT. 

An interesting feature of equations 22 and 23 is that 
< S z > diverges as the square of the logarithm of the 
cutoff, signaling that the Renormalization Group analysis 
does not work at that point. 



V. CONCLUSION 

The purpose of this paper was to shed some new light 
on the recent controversy about the existence of a spin 
liquid phase in the phase diagram of the J% — J 2 and 
•h — Ji — J3 models. We have tried to do a thorough 
analysis of the problem using different techniques, com- 
paring them and extracting consequences about the va- 
lidity of each one. We have obtained a very consistent 
qualitative picture of how the phase diagram looks like. 
For the J\ — J 2 model, there seems to be a first or- 
der phase transition, due to the fact that the stability of 
the Neel state is enhanced by quantum fluctuations. At 
the same time there probably exists a very small region 
where the magnet is disordered for spin i. Its frustrated 
point, though, is a Lifshitz point and therefore is not 
representative of the phase diagram of the J\ — J 2 — J3 
model. In the J\ — J2 — J3 model we have found or- 
der from disorder, too. But now there is a finite region 
where the system is disordered for any value of the spin, 
because the stability of the helicoidal phase is reduced 
as much as the one of the Neel phase is enhanced, and 
there is no first order phase transition. For spin ^ the 
spin liquid region is shifted from the Classical Critical 
line Ji — 2 J 2 — 4 J 3 = (figure 1) to considerably bigger 
values of J 2 and J3. We find that the character of the 
spin disordered phase seems to be different for the weakly 
and for the strongly frustrated cases. This fact could 
have any repercussion in the role of instanton tunneling 
between different topology sectors, and the stabilization 
of spin-Peierls order close to the critical line. Finally, we 
would like to mention that our results agree with numer- 
ical diagonalization studies by Dagotto et al 1||l3l both 
for the J\ — J2 — J3 and the J\ — J2 models but are in 
contradiction with the recent work of Schulz and Ziman 
[14 1 and with the results obtained from series expansions 
[12 1 for the J\ — J2 model. 
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APPENDIX A: SPIN WAVE THEORY FOR AN 
HELICAL MAGNET 

For a general Heisenberg model 

n = J2 J ij^i-s j (Ai) 

a classical helical configuration Si — S(usm(Q ■ Ri) + 
v cos(Q • Ri)), with u ■ v = 0, has an energy E c i — 
N S 2 J(Q). This helical state lies in the plane normal 
to the "twist" vector n = u x v. The classical ground 
state is the helix with pitch Qo.cZ obtained by minimiz- 
ing E c i, i.e.: 9^J(<5o,c/) = 0. To perform a spin wave 
expansion about this classical state we rotate the local 
quantization axis through an angle 9i = Q ■ Ri about the 
twist vector h, transforming to a twisted reference frame 
in which the ordered state is ferromagnetically aligned 

IB 

C. __. pAiX q. _ Q RiflX q. 

with this prescription, we obtain the following Hamilto- 
nian: 



H ' = Yj J v {"f?A ■ + 7^(2 Si-nS r h-Si- Sj) - sin(Q • R^fi ■ (S t x S,)} (A2) 



ij 



where 7y = h (1 ± cos(Q ■ Rij))- We choose n — z to 
carry out the calculations of this section. 

The usual way to deal with these magnetic models 
is to decompose the spins into Holstein-Primakoff (HP) 
bosons |58],[39| and treat the resulting Hamiltonian as a 
dilute gas of bosons Q . We use this representation to 
calculate LSWT properties of the different spin configu- 
rations we are dealing with. The ground state energy of 
the helicoidal Q = (q, q) state is 



E °- S{S+1) jQ + ^E(k) (A3) 



N 2 vo 2 

k 



E i k ) = \j (JQo - J k)(JQ ~ ^(Jk+Q + Jk-Q )) = \J (Jq ~ Jk)(JQ ~ J±) 



where again Qq is determined by minimizing Eq. We 
can expand Qq in powers of 1/S about its classical value 

(<7o : d = " 2j 2 +4,/ 3 for tlie state Qa = go))- The result 
is 

1/2 



*-4"- A 4-5«U?*'*(^|) +0( i» (A4) 



The spin stiffness is defined as the second derivative 
of the free energy with respect to a uniform twist to the 
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direction of the order parameter |4l| . This is equivalent 
to the response function to a gauge field Ai = 5Q ■ Rik 
applied to the system, where k is the axis of the twist 
and SQ, its wavevector. Then, 

d 2 F 

(A5) 



dAfdA"; 



For an helical magnet, with n || z, we have p^ = p v u 7^ p\ { 
and p xx = Py y . For the purpose of this paper, it will be 
enough to calculate pf,. We obtain: 

z _ S 2 diiJQ acl 
Pa - 2 

+ 2 {9uJ qo , cl - (sin(«o,d) + 8 J 3 sin(2Q ci ))SA^ 

1 \ - (2Jq 0iC! - Jk - J±)d li JQ acl - (Jq , c! - Jk)duJ± 1 

On the other hand, for a bipartite lattice, and when 
the ground state is Neel ordered (n, 7r), it proves more ad- 
vantageous to use a Dyson-Maleev (DM) decomposition 
of the spins |38| , because then, the Hamiltonian doesn't 
have terms higher than fourth order and because cor- 
rections to mean field results are less important than the 
ones obtained with HP bosons. The results obtained with 
this method are identical to those of HP at the mean field 
level as we will show p0[ . 

We want to calculate the staggered magnetization for 
the Neel state of the two models to next to LSWT or- 
der in 1/S using DM to check the convergence of the 
series. The Hamiltonian is given by equation A2, with 
Q = (tt,tt). We perform a Fourier transform after intro- 
ducing DM bosons and rewrite it as: 

H = SzJ^ i A (k) {a+a k + b+b k ) + J llk (a k b_ k + a+bt k )} 

k 

- ^ 5(1 + 2 - 3 - 4) {2 7 1+2 a+a_ 2 ^ 3 &4 + li<4b+ 6± 3 6 4 + l4flta± 2 

1234 



" h ]T 6(1 + 2 - 3 - 4) {(r^ - T A 2 _ 4 )a+a+a+a 4 + (T* - T s 2 _ 4 )b+ b+ b 3 b 



N 



where 



A(k) = J 1 -J s (l-T s k ) 

Ik = 2 (cos(k x a) + cos(k y a)) 
T k = cos(k x a) cos(kya) 

rl = 72fe 

and z is the coordination number. 

We introduce a Bogoliubov transformation, = 
ua^ + vj3- k and b_ k = vat 4- uf3- k , and normal-order 
all terms of the Hamiltonian with respect to the vac- 
uum of a and j3, neglect the contribution coming from 
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normal-ordered terms involving four Bogoliubov bosons 
and end up with a quadratic diagonal Hamiltonian. The 
procedure is standard E]. We obtain a self-consistent 
equation for the Bogoliubov amplitude Vk, by requiring 
that the terms proportional to a + (3 + and a/3 are zero: 

Vik) = - S ^^\ F m-l\ m (A8) 



V2 I £{k) 
2 

TV 



Fx(fc) = A{k)s - 1 {v 2 (fc')(i - Mi - rl - it + r L*0) + iwUtfW)} 



k' 



F 2 (k) = lk S + bfcV(fc') 2 + 7fc+ fc ^(fc')V(fc')} 



TV 

k' 



£(k) = yjF?(k)-F*(k) 



here £ (fc) is the energy of the spin waves and we have 
used units where J\ = 1. This expressions can be ex- 
panded in powers of 1/S. Doing so we obtain the follow- 
ing relation for the staggered magnetization in terms of 
LSWT functions: 

c i 1 1 y A W J * y ^ 2 (i - rf) v A(fc')it - 7fc 2 , fAql 

" " 2 N ^ E(k) N 2 S 4^ £^ 3 (fc) V 



where 



E{k) is the energy of t he spin waves in LSWT. We note 



here that equation A9 agrees with the expression for the 
staggered magnetization obtained in reference 2. 

We can also compute £(k) to next to LSWT order in 
1/S in the long wavelength limit in terms of the spin 
wave velocity 

c a (5)^8Ji5o|poU + |(l + ai-a4))-^(o2-ai)J> (•/.■; = i>i (AID) 
-SJiSajpo U + |(l + ai-a 4 )J -^(a 3 -ai) J- (J 2 = D) 



"1 = > , TT7TT «2 = 2^ ■ 



^£(fc) - ^ £7(fc) 



A{k)Tl ^ A(k) 

a3 = ?^T ai ^W) 



APPENDIX B: SELFCONSISTENT EQUATIONS 
FOR SCHWINGER BOSONS MEAN FIELD 
THEORY 

Some years ago, Arovas and Auerbach jl3| introduced 
a SU(N) functional integral expansion as a useful tool 
to deal with unfrustrated magnets which are disordered 
due to quantum or thermal fluctuations. It was based on 
a SU{N) generalization of the spin decomposition 
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§i=~bfffbi bfbi = 2S bt = (6+,6+) 

originally due to Schwinger f3S|] . We will use here the 
generalization to frustrated magnets due to Coleman and 

coworkers H@ for S^ 2 )- 

We begin by rotating the local quantization axis an 

angle 6% — Q ■ Ri so that in the twisted reference frame 
the spins are ferromagnetically aligned, and introduce the 
triplet Cooper and singlet particle-hole even-parity fields 

B% = /bVb; D+ = b+b, 

with which, the action reads (choosing h = z): 

S = So dT { h+drh ~ 2SNX ~ NS2 2 JQ + XD + \ DJ+D - \b z+ J-B z ) i Bl ) 

J kk> = 2 ± ^k-k'+Q + Jk-k'-o) 

where summation over k,k' is understood, and \{t) 
is a Lagrange multiplier which enforces the constraint at 
all sites and times, and which can be regarded as another 
field. 

We perform a Hubbard-Stratonovich transformation 
from Dk, Bk to the bond fields hk, A& and make a saddle- 
point approximation of the resulting action in those bond 
fields. At the saddle-point, they are time independent, 
spatially uniform and real (and so is the field A). After 
diagonalizing the resulting mean field Hamiltonian with 
a Bogoliubov transformation bk a = wjker + v l-k-a wc 
find the following Free Energy 

£ = -(25 + 1)A - ^ - ^ J+T 1 /^ + A, J"- 1 A V + - hi ( 2 sinh ^ ) , B2 , 



e(fc) - V(Mfc) + A) 2 -A 2 (fc) 

When the ground state is ordered, we have to allow for 
a one particle condensate contribution to the mean field 
equations glUHl. Besides, the gap in the spin wave 
spectrum closes -we have a Goldstone mode-, so that we 
have to adjust A to fulfill this requirement. The saddle 
point of F at T = is given by the following equations: 



£ = _(2S + 1)A - ^ - £ + E P* ( B3 ) 



at = S*Sko + 

Pk = S*5kQ + 



hk + A 
2e(fc) 
A fc 

2e(fc) 



e(fc) = v/(^-A) 2 -A2(fc) 

In general, the set of values of the pitch Q for the 
ground state are obtained by minimizing F with respect 
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to Q. We can obtain F as a series in 1/S by iterating the 
mean field equations. We find 



e(k) = S^Jq - J k )(J Q - -{J k+Q + J fc _ Q )) (B4) 
F^ + S J Q +J2e(k) = F S 2 + F 1 S 

k 

The expression for the spin wave energy coincides with 
that given by LSWT while F± is twice the value given by 
LSWT (as already stated by Arovas and Auerbach B3] 
for the unfrustrated case, SBMFT overcounts the number 
of degrees of freedom) . 

SBMFT has two major problems in dealing with he- 
licoidal states: (1) Symmetry considerations require the 
existence of Goldstone bosons at (0,0) and (<?, q). We 
can adjust A so that we have one of them, say (0, 0). But 
then, we don't have the other one: \hQ + Aq — ho\ =^ |Aq| 
(for Neel and collinear states we do get both zero modes) . 
This leads to unphysical results for the stability diagram 
of the helicoidal phase, and in particular, to a first order 
phase transition. (2) We have already shown that F\ the 
contribution of order S to the free energy is twice the 
correct value. This means that the pitch Q obtained in 
the minimization is not the correct one. 



APPENDIX C: DERIVATION OF THE SCALING 
EQUATIONS FOR THE LONG WAVELENGTH 
ACTION 

In this appendix we will use the momentum-shell 
method to derive recursion relations for the S-matrix of 
the hydrodynamic action of the J\ — J2 — J3 model, equa- 
tion This is a generalized 0(3)/f/(l) = S 2 model. 
The unit vector n — (it, er); tt is a n — 1 = 2-component 
vector, and a is the direction of quantization. The S- 
matrix is given by 

S(-oo,oo) = I VhY[6{a 2 {x,t) +n 2 (x,t) - l)e lS 

J x,t 

(CI) 

It is convenient to define the dimensionless variables r = 
r/a and t = ct/a so that the action is 

-1 />oo />L 

S = — / dld 2 f {(^n) 2 -((^) 2 ) + (%n) 2 )) (C2) 

+ 12|) (( ^" )2 + W"*> a > - ^-f^)(^n) + 2 -^f] 

We integrate a(x, t) now, and generate an action which 
is nonlinear in the fields 7r(x, t) . The nonlinear terms 
constitute the interaction p4| . As we are interested in 
one loop corrections to the charge g, we retain only the 
quartic terms. We divide the action now into a free part 
and the interaction 
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1 f dwd 2 k 

5 i = ( II / • *( 2 ) -( 3 ) ' + 2 + 3 + 4) 

x |w 2 W4 -k 2 -k 4 - ^ + j^^(k xl k x2 k x3 k x4 + k 2 x2 k 2 xA + 2k xl k x2 kl 4 + (x <-> 
-j^{k x ik x2 k yi k yA + k 2 x2 k 2 y4 + 2k xl k x2 k y4 )^ + ^ ^^^ (fc, tu) • 7r(-fc, - 



where /i = " gg and A t is an ultraviolet cutoff for the 
time direction. The free action defines an unperturbed 
propagator Qo, which we use to compute the different 
terms which arise in the loop expansion. We compute 
the renormalization of the free term to one loop order, 
by expanding e Sl , contracting fields and exponentiant- 
ing again. The measure cancels the ultraviolet divergent 
terms arising in the one loop expansion as it does for the 
classical case. In particular, the regulator A t doesn't ap- 
pear in the renormalized theory. The procedure is very 
well explained in reference 30, so we refer the reader to 
it for more details. We obtain 



x j W 2 (l + g l) fc 2 (l + .90/(1 - )) + 12^5) ^ + W + ^) 

- -^klk 2 y (l + g I)-goh(l + goI)} 



where I is the following loop integral 



( ( dw f 1 d 2 k 



2tt J e -> (2tt) 2 w 2 k 2 + 12 p(g) (^ 2 + k y ) ~^s)k 2 k 2 
v / 6p(S)l r /2 d9 



Ms) i r 

* 2 Jo 



\J 2Ap{S) - 2a i + (cti + 6cr 2 ) sin 2 6 



It diverges at the tricritical point, showing that the 
scaling analysis can not be applied to dilucidate what 
happens in the J\ — J 2 model close to the tricritical 
point. 

After performing the rescaling 

c > 

k = e~ l k' w = —e~ l w' 
c 

TT = (n' ti = Ch 

we find the parameters of the renormalized action in 
terms of the old ones: 

(ghy = e 2l (gh)(l + gl^^j (C5) 
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(-)' = -( 2 e~ 51 (l + gl(l 3ff l^ 



<T V ' 8 p(S) 

c[l-gl- 



8p(S) 
<Ti = <n(l - 21 - gl) 



,' (S ) = ,( S )(i- 9 /(i + ^|f)) 

from which one obtains easily the scaling equations of the 
model after iterating. 
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FIG. 1 . Classical Phase diagram for the Ji — J2 — J3 model. 
Thick solid lines indicate a discontinuous transition; dashed 
lines indicate a continuous one. The thin lines are the results 
of LSWT for the boundaries of the different ordered states for 
spin 1/2. 

FIG. 2. Phase diagram for the Ji — J2 — J3 model. Dashed 
lines are LSWT; solid, SBMFT. The SWT p Neel vanishes 
along the dotted line. 

FIG. 3. (a) LSWT correction to go jC ; as defined in equa- 
tion 8. (b) Classical spin stiffness (solid lines) and LSWT 
corrections to it (dashed lines). 
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FIG. 4. (a) Flow diagram for the QNLctM in 2+1 dimen- 
sions. The dashed line separates the Renormalized Classical 
region from the Quantum Disordered one. (b) Flow diagram 
for the Ji — J2 — J3 model, linking CHN and Ioffe and Larkin 
analysis. We plot in the y-axis p(S), the spin stiffness renor- 
malized by short-wavelength quantum fluctuations. Rotating 
the figure by 90° we obtain the phase diagram of SWT. 

FIG. 5. (a) Phase diagram for the Ji — J2 model. Dashed 
lines are LSWT results, while solid lines are those given by 
SBMFT. (b) same as before in a logarithmic scale. The cross 
is the result from numerical diagonalization of small clusters 
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